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The population dynamics of a trapped Bose-Einstein condensate, subject to the action of an oscil- 

. . , latory field, is studied. This field produces a modulation of the trapping potential with the frequency 

• close to the transition frequency between the ground state and an excited energy level. Unusual 
I critical effects are found exhibiting sharp qualitative changes in the population dynamics. An ef- 

• fective averaged system is constructed explicitly illustrating the occurrence of critical phenomena. 
' The related critical indices are calculated. 
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Bose-Einstein condensates of trapped atoms present an interesting example of essentially nonlinear statistical sys- 
' terns with rich properties that are being intensively studied both experimentally and theoretically [1,2]. Although 
] experiments deal with rather dilute Bose gases, their interaction, nevertheless, cannot be considered as a small per- 

■ turbation. This is because in a system of condensed atoms coherence develops, resulting in an effective interaction 
' proportional to the number of particles. The fact that Bose-Einstein condensation and coherence appear simultane- 
. ously can be naturally understood remembering that both these phenomena arise when the thermal wavelength of 
' atoms exceeds the mean interatomic distance. All the measurements done so far show the evidence for coherence of 

Bose-Einstein condensates [1,2]. Essential nonlinearity of the latter makes their behaviour quite different from that 

■ of ensembles of weakly interacting gases. 
In the present article we show that the nonlinearity in the Hamiltonian describing Bose condensates results in 

very unusual dynamical effects resembling critical phenomena in statistical systems. We consider the dynamics of 
' population levels of a Bose condensate subject to the action of a resonant field whose frequency is close to the transition 
. frequency between the ground-state and a chosen excited level. Varying the characteristic parameters, such as the 
(~| ■ transition amplitude and detuning, we show that there exists a bifurcation line, dividing the space of these parameters 
O I onto regions, where the dynamical behaviour of the system is dramatically different. This bifurcation line is analogous 
O ' to the critical line of a statistical system experiencing critical phenomena. To substantiate this analogy, we construct 
^ , an effective system explicitly demonstrating critical behaviour along the critical line coinciding with the bifurcation 
line found for the evolution equations. For the effective system, it is possible to define the quantities analogous to the 
order parameter, specific heat, and susceptibility and even to calculate their critical indices, 
i Let us consider a Bose system, when all atoms are condensed and the system can be described by the nonlinear 
■ " ' Schrodinger equation, often called the Gross-Pitacvskii equation [1,2]. Assume that, in addition to a stationary 
trapping potential, there is a time-dependent potential, which we call the resonance field and whose meaning will be 
specified below. Thus, we consider the equation 
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with the nonlinear Hamiltonian in which [/ is a trapping potential; A = ATrh^asN/mo; Ug is a scattering length; mg 
is mass; and N is the number of atoms. The resonance field has the form Vres = V{f' ) cosojt. This can be treated as 
the modulation of the trapping potential. 

Tuning the frequency to of the resonant field close to the transition frequency cumn = {Em — En)/h between the 
energy levels defined by the stationary problem H{ipn)(pn = E„ipn, we may realize the corresponding interlevel 
transitions [3]. Since the nonlinear Schrodinger equation is an exact equation for the wave function of a coherent 
state [4], the stationary solutions (pn can be called coherent modes. These should not be confused with collective 
excitations defined as small deviations from the ground-state and described by a linearized equation [1,2], although 
resonant transitions between different branches of collective excitations are also possible [5]. But the equation for 
coherent modes is principally nonlinear. 
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Suppose that at the initial time t = the Bose gas was condensed, with ah atoms being in the ground state 
(/?(r, 0) = ifoij). After this, the resonance field is switched on, with a frequency lu being in quasiresonance with the 
transition frequency ojjo between the ground-state level and a chosen energy level j. The quasiresonance condition 
is |Au;| ^ iVjo, where the detuning Aw = u — ojjo- The solution of Eq. (1) can be presented as a sum ip{f,t) = 
J2„ Cn{t) fn{f,t) over coherent modes (/?„(r, t) = (pn{r) exp{—iEnt/h). Substituting this expansion in Eq. (1), one 
gets a system of equations for the coefficients c„ (t) . These equations can be simplified employing the quasiresonance 
condition. To this end, it is convenient to introduce the notation for the fractional level population ni{t) = |ci(t)p, 
and for the interaction intensity a and the transition amplitude f3, respectively. 



for the dynamics of the groimd-state coefficient co{t) and of the coefficient Cj(t) related to a chosen excited mode j. 
According to the assumption that at the initial time all atoms are condensed in the ground state (po{f), the initial 
conditions to Eqs. (3) are Co(0) = 1 and Cj(0) = 0. The derivation of Eqs. (3) has been explained in detail in Ref. 
[3]. However, it is worth emphasizing here a couple of important points. Looking for a solution of Eq. (1) in the 
form of the sum Cntpn does not require that the modes (pn compose a complete set. Recall that the completeness 
of a basis is the possibility of expanding over it an arbitrary function from the considered space. Such an excessively 
restrictive property is of no need for us. What we need is just to present the sole function </? as a sum over modes, 
with the coefficients c„ to be defined from Eq. (1). In the theory of differential equations, this way of looking for 
a solution is called the method of substitution or the method of variation of parameters. This method is often used 
for solving nonlinear differential equations, e.g. for nonlinear optical waveguide equations [6,7] and for the nonlinear 
Schrodinger equation [8,9]. An exact orthogonality of the nonlinear modes is also not necessary, but, in general, it is 
sufficient that the modes tpn be approximately orthogonal in the sense of the smallness of the scalar product {(pm, pn) 
for n. In the case of the nonlinear Schrodinger equation, taking for coherent modes the variational wavefunctions 
of Ref. [3], it is easy to check, by direct calculations, that the maximal values of \{prm'-Pn)\ for different m ^ n are 
of order 0.1. Moreover, even this approximate orthogonality is not compulsory in the frame of the quasiresonance 
picture employed in Ref. [3]. For deriving Eqs. (3), it is sufficient that c„(i) be considered as a slow function as 
compared to exp(— i£J„f/7i). Then, integrating Eq. (1), one can make use of the exact orthonormality on average of 
functions pn in the sence of the equality limr^oo 7 Jq [J P*n{^7 t)fn{r, t)dr'\ dt = Smn- Calculations of Ref. [3] show 
that Cn{t) can, really, be treated as slow, compared to exjp(—iEnt/h) since the variation rate of c„ is of order a which 
is an order smaller than i?„, that is, \dc„/dt\ <C En- 

We solved the system of nonlinear equations (3) numerically, carefully analyzing the behaviour of the solutions 
for different parameters. This behaviour turned out to be surprisingly rich exhibiting unexpected critical effects. 
Analyzing Eqs. (3), it is convenient to make there a scaling, measuring time in units of and introducing the 
dimensionless parameters b = \P\/a and S = Au/a. The dimensionless detuning is assumed to be always small, 
5 <C 1. When the dimensionless transition amplitude b is also small, the fractional populations oscillate according to 
the sine-squared law. When b increases, the amplitude of oscillations also increases. The overall behaviour continues 
to be normal unless we reach a critical value be, when the dynamics of the system changes drastically. More generally, 
there exists a critical line connecting the parameters b and 6, so that b^ = 0.5 — 5. In the vicinity of this line, the 
system dynamics experiences sharp changes, when the parameters are varied just a little. We illustrate this in Figs. 
1-4, slightly varying the detuning and keeping b — 0.4999, so that we are close to the critical line. In Fig. 1, the 
detuning is zero, and the oscillations of the fractional populations are yet normal. Shifting the detuning to S = 0.0001 
transforms the picture to that in Fig. 2, where the top of nj{t) and the bottom of no{t) become flat, while the 
oscillation period is approximately doubled. Changing the detuning to 6 = 0.0001001 results in Fig. 3, where the 
period is again doubled and there appear the upward cusps of nj{t) and the downward cusps of no{t). Increasing 
further the detuning to 5 = 0.00011 squeezes the oscillation period twice, as is shown in Fig. 4. The same phenomena 
occur when we cross the critical line 6 + (5 = 0.5 at other values of parameters. 

This unusual behaviour of the fractional populations is certainly due to the nonlinearity of the evolution equations 
(3). Systems of nonlinear differential equations, as is known, can possess qualitatively different solutions for parameters 
differing by infinitesimally small values. The transfer from one type of solutions to another type, in the theory 
of dynamical systems, is termed bifurcation. At a bifurcation point, dynamical system is structurally unstable. 
Bifurcations in dynamical systems are somewhat analogous to phase transitions and critical phenomena in equilibrium 
statistical systems. To elucidate this analogy for the present case, we have to consider the time-averaged behaviour 





Then, in the quasiresonance approximation [3], we have the system of equations 
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of the system, which can be done as follows. First, we need to define an effective Hamiltonian generating the 
evolution equations (3). To this end, we notice that the latter equations can be presented in the Hamiltonian form 
idco/dt = dHeff/dcQ, idcj/dt = dHeff/dCj, with the effective Hamiltonian 

H^ff = a nonj + ^ (/? e^'^'^^^cj + p* e-^^"*c*co) . (4) 
Then, using the averaging method, we solve the evolution equations (3) finding 

expj- ^(a- Aw) t| , 

Cj = — i — sin — exp i — -{a + Aw) t>, = [a{no — rij) — Aw] + |/3|^ . (5) 

The usage of the averaging method here is accomplished in the same way as it has been done for solving nonlinear 
evolution equations describing superradiant spin relaxation in magnets [10,11], or nonadiabatic dynamics of atoms in 
magnetic traps [12,13]. Technical details can be found in the latter references. Employing the found solutions (5), with 
the normalization ng + rij = 1, we obtain for the effective Hamiltonian He// = a + rij Aw. Averaging the fractional 
populations over the explicitly entering time, we get the average population Uj = b/2e, where, as earlier, b = \(3\^/a 
and e = ft/ais the dimensionless average frequency, with il given in Eqs. (5), where n is replaced by rij. This average 
frequency is defined by the equation e^{e'^ — iP') = (e^ — 6^ — e'^S)'^. Substituting the average population nj in the 
effective Hamiltonian, we obtain the average effective energy -Eg// = (Q;6^/2e^)(6^/2e^ + S). An order parameter for 
the averaged system can be defined as the difference of the average populations, r] = no — nj = l — h'^/e^. The capacity 
of the system to store the energy pumped in by the resonant field is characterized by the derivative = dE^j j / d\P\ 
of the average energy with respect to |/3|. And the absolute value of the derivative of the order parameter with respect 
to 5 defines a kind of susceptibility xs = l^^/^^l- Note that plays the role of heat capacity, because of which it 
can be called the pumping capacity. 

The solution for e shows that the average frequency changes from e = l — 2(5at5=:0toe = \/0.5 — 8 at 
b = be = 0.5 — 5. After this, the frequency diminishes by a jump to e = 6 for b > be- This implies that at the 
point be the average oscillation period increases by a jump almost twice. The order parameter diminishes from 
r] — 1, when there is no resonance modulation, that is 6 = 0, to rjc = 0.5 + 6 a.t b = be, after which the order 
parameter becomes approximately zero for b > be- The pmnping capacity grows from zero at & = 0, becoming 
divergent at b = be- The same concerns the susceptibility. All this suggests that be is a critical point and the 
relation be + 5 = 0.5 defines the critical line, which is in agreement with our numerical solution of the evolution 
equations (3). The asymptotic behaviour of the characteristic quantities in the vicinity of the critical line defines the 
corresponding critical indices. Thus, for the small relative deviation t = \b — be\/be from the critical point be, 
we have rj — r]e — ^ {I - 2(5)t^/^, C/j ~ t~^/^, xj ~ r"-^/^. Hence, the related critical indices are 1/2 
for all characteristics. It is interesting that the critical indices of rj, Cp, and xs satisfy the known scaling relation: 
ind(Cfj) + 2 ind{r]) + ind{xs) = 2, where ind is the evident abbreviation for index. 

If, after exciting a coherent mode, one switches off the pumping, then the following dynamics is to be considered 
separately, taking into account stability conditions as well as the lifetime of the condensate as such, which is defined 
by depolarizing collisions [14]. Finally, let us analyse the optimal conditions when the considered effects could be 
achieved experimentally. To reach the bifurcation line one needs to invoke the transition amplitude /3 « a/2. At the 
same time, the interaction parameter a, at least in the strong-coupling limit [3], can be of the order of the transition 
frequency w. When the magnitude of the transition amplitude [5 becomes close to the transition frequency w, the 
effect of power broadening is expected. This effect is well known in optics [15], where the role of /3 is played by 
the Rabi frequency. In the regime of power broadening, one would expect the ground state to be coupled to more 
than one excited mode, so that the quasi-resonant two-mode picture could become not a good approximation. This 
picture remains a good approximation provided that the probability of nonresonant excitation of neighboring levels is 
small. The probability of an induced transition between levels i and j, due to a monochromatic field with frequency 
w, can be estimated [15] as « i3fj/2[{uj — uJij)^ + fif^], where (iij and w^ are the related transition amplitude and 
frequency, respectively. To estimate , we may use the calculations of Ref. [3] , slightly modifying them to take into 
account the cylindrical symmetry of a trap with radial and axial frequencies w^ and ujz, respectively. Dealing with 
a cylindrical trap gives an additional explicit parameter v = Wz/uJr, hence, more possibilities for varying conditions. 
In the case of nonresonant excitation, when w — Wy ^ w, the magnitude of essentially depends on the value of 
the transition amplitude /Jjj. The maximal value of the latter, occurring when approaching the bifurcation line, is 
close to a/2. The value of a depends on the parameter gv, with g = AiragN/lr, where is a scattering length, N is 
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the number of atoms, and Ir = y/h/motOr. If gi^ is not large, then a <C t^. Thence the maximal Pij « a/2 <C co and 
Pij <C 1. Therefore the nonresonant excitation can be neglected, and the two- level picture is a good approximation. 
When gi/ ^ 1, then a ^ lu, so that the maximal (3ij ~ i^/2, from where Pij ~ 0.1. This tells that, although power 
broadening can influence quantitative results, the two- level picture still serves as a reasonable first approximation. To 
reduce the influence of power broadening, one can proceed as follows. In the case of an almost spherical trap, with 
^' « 1, the number of atoms is to be restricted in order not to make g too high. When the number of atoms is large, 
so that g ^ 1, this can be compensated by taking a long cigar-shape trap, with <C 1, which would result in gi/ ~ 1. 
In the present-day experiments there is such a wide variety of possibilities for changing characteristic parameters [1,2] 
by choosing an appropriate shape of a trap, by changing the number or the type of trapped atoms, by varying the 
scattering length using Feshbach resonances, and so on, that it looks rather feasible to observe the critical phenomena 
described. 
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Figure captions 



Fig. 1. The time dependence of the fractional populations no{t) and nj{t) for b = 0.4999 and 6 = 0. Here and 
in all following pictures the dashed line corresponds to the ground-state population no{t) and the solid line, to the 
excited-level population rij (t) . 

Fig. 2. Flattening of the fractional populations, with their oscillation period being doubled, at 6 = 0.4999 and 
6 = 0.0001. 

Fig. 3. The appearance of the upward cusps of nj{t) and of the downward cusps of no{t) for b = 0.4999 and 
6 = 0.0001001. 

Fig. 4. Fractional populations versus time for b = 0.4999 and d = 0.00011. 
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